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Abstract: We derive the quadratic action of the fluctuations around the classical world sheet 
underlying the jet quenching from AdS/CFT. After obtaining the correspondence partition 
function, the expansion of the jet quenching parameter of A/" = 4 super symmetric Yang- 
Mills theory is carried out to the sub-leading term in the large 't Hooft coupling A at a 
nonzero temperature. The strong coupling corresponds to the semi-classical expansion of 
the string-sigma model, the gravity dual of the Wilson loop operator, with the sub-leading 
term expressed in terms of functional determinants of fluctuations. The contribution of these 
determinants are evaluated numerically. We find the jet quenching parameter is reduced due 
to world sheet fluctuations by a factor (1 — 1.97A~^/^). 
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1 Introduction 

One of the main purposes of the heavy-ion colUsion experiments is to explore the QCD 
phase diagram and the properties of new state of matter created through collisions at high 
energy density. The experimental program at the Relativistic Heavy Ion Collider (RHIC) in 
Brookhaven National Laboratory has led to strong evidences that this new state of matter is 
in the deconfinement phase of QCD, but behaves very differently from a free gas of partons. 
The name "strongly coupled quark-gluon plasma (SQGP)" was given to this matter because 
of several striking experimental discoveries. Firstly, the measured elliptic flow data can be 
well described by hydro dynamical models only with very low shear viscosity, suggesting that 
SQGP resembles a nearly perfect liquid. Secondly, due to the interaction with the medium, 
high energy partons traversing the medium are strongly quenched. This phenomenon is 
usually characterized by the so-called jet quenching parameter (or transport coefficient) q, 
which describes the average transverse momentum square transferred from the traversing 
parton, per unit mean free path [1, 2]. There are still some model-dependent uncertainties 
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on the extracted values data of the transport coefficient from heavy- ion colhsions, ranging 
from 1 — )• 25GeV'^ / fm [3, 4] , which is considerably larger than those from perturbative QCD 
estimation [5, 6], higher values demand additional non-perturbative mechanisms. Therefore, 
it is of great importance to get further information on the possible values of q in the strong 
coupling limit. The conventional nonperturbative first principle tool namely lattice QCD, 
cannot be applied for this purpose, which requires real-time dynamics. 

A prominent implication of the AdS/CFT duality [7-10] is the correspondence between 
the type IIB superstring theory formulated on AdSs x and A/" = 4 supersymmetric Yang- 
Mills theory (SYM) in four dimensions. In particular, the supergravity limit of the string 
theory corresponds to the leading behavior of SYM at large Nc and large 't Hooft coupling A = 
Sym^c = ^ with L the AdS radius and a' the reciprocal of the string tension. This relation 
thereby provides a powerful tool to explore the strongly coupled QGP created in RHIC in 
spite of its underlying dynamics, QCD, is different from A/" = 4 SYM. It is expected that some 
of the properties of the latter is universal for all strongly interacting system. Computationally, 
AdS/CFT correspondence translates the quantum field theory of supersymmetric Yang-Mills 
theories at strong 't Hooft coupling and large number of colors Nc into a classical field theory 
in a gravitational background. The AdS/CFT correspondence has been successfully applied 
into the RHIC physics, including the thermodynamical the transport properties of sQGP. One 
of the celebrated finding is the universal value of the ratio between the shear viscosity and the 
entropy density, rj/s = for quantum field theories admitting a holographic description. It 
was further conjectured that rj/s > ^ for all matters. See [11] for a review on the subject 
and the references therein. The RHIC data of the viscosity ratio is compatible with this lower 
bound. In addition, the entropy density of a = 4 SYM plasma normalized by its Stefan- 
Boltzmann limit equals to 3/4 [12], which is close to the corresponding quantity extracted 
from lattice QCD. 

Motivated by these similarities, Liu, Rajagopal and Wiedemann (LRW) computed the 
jet-quenching parameter of a AA = 4 SYM plasma with the aid of AdS/CFT correspondence. 
The Eikonal approximation relates the jet quenching parameter with the expectation value 
of an adjoint Wilson loop 1^^[C] with C a rectangular contour of size L x L_, where the sides 
with length L_ run along the light-cone and the limit L_ — t- oo is taken in the end [13]. Under 
the dipole approximation, which is valid for small transverse separation L, the jet-quenching 
parameter q defined in Ref.[l] is extracted from the asymptotic expression for TL << 1 

< W^[C] Woexp[-^qL^L^] (1.1) 

with T the temperature and the prefactor Wq independent of L. The correspondence principle 
relates < VF"^[C] > to the quantum effective action of a type IIB superstring in AdS^ x 
with a black hole. In the strong coupling limit, the action is proportional to the minimum 
area of the string world sheet in the AdS^ bulk spanned by the C at the boundary. They find 
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that [14] 

Interestingly the magnitude of q^'yAi turns out to be closer to the value extracted from RHIC 
data [15, 16] than pQCD result for the typical value of the 't Hooft coupling, A ~ Gvr, of QCD. 
This proposal has attracted lots of interest. Alternately, the energy loss and jet quenching 
problem was also studied from a drag force on the brane by several authors [17-21]. Their 
extracted jet quenching parameter is also proportional to ^/XT^ but with a smaller coefficient 
(by 20 percent). 

The LRW formula is strictly valid only when the 't Hooft coupling constant goes to infinity 
and large Nc limit. However, with a finite value of A and the Nc = 3, an understanding of 
how these computations are affected by finite A corrections may be essential for more precise 
theoretical predictions. There are two sources of contributions to finite A corrections: 1) 
The string theory corrections of the AdS-Schwarzschild background where the world sheet 
underlying the jet quenching Wilson loop is embedded. 2) The fiuctuation of the world sheet 
itself. The former has been calculated in [22] and it amounts to multiply (fgyM ^ factor 
1 — 1.765A^^/^. The parallel correction within the framework of the drag force has also 
been figured out in [18], and has been found to enhance the jet quenching parameter. The 
latter contribution, that from the fluctuations, however, is of the order 0(A~^/^) and thereby 
dominates the former. But it is much harder to evaluate. It is the purpose of the present 
paper to determine this contribution and we obtain the corrected jet-quenching parameter 

qsYM{X) = 4yAf [1 + + 0(^-1)] (13) 

with the value of k 

K~-1.97. (1.4) 

The paper is organized as follows. In the next section, we will derive the partition 
function with the world fluctuations. We will present our analytical study and numerical 
results in section 3 and section 4 respectively. The section 5 concludes the paper along with 
some discussions of the result and some open issues. Some calculation technique details are 
presented in the appendixes. 



2 The one loop effective action 

It follows from the AdS/CFT correspondence that the expectation value of a Wilson loop 
of the N = A SYM in the fundamental representation of the gauge group is related to the 
path integral of a type IIB superstring in a 10-dimensional spacetime that is asymptotically 
AdS^ X i. e. 

< W[C\ >= const. [ [dX][de]e*^[^'^l = e^^-^'P], (2.1) 
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where X's and ^'s are the bosonic coordinates and fermionic coordinates, and the superstring 
action 

s[x,e] = SB[x] + SF[e] (2.2) 

with 5b [X] and Sf[0] the bosonic and fermionic parts [23, 24]. The semiclassical approxima- 
tion of the effective action ^cff. reads 



5eff.[C] = ^[So[C] + + O ] , (2.3) 



where \/A<S'o[C] = 5'[X,0] with X the classical solution and Si[C] stands for the one-loop 
correction. With the light-like antiparallel lines, the jet-quenching parameter comes from the 
second power of the Taylor expansion of Scs. in 
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(27r)2TL , 

(2.4, 

with L the distance between the two lines and T the temperature. 

To figure out >S'i[C], we need to expand S[X, 9] to the quadratic order around the classical 
solution X and ^ = 0. To the same order of A, the 10-dimensional metric takes the form of 
Schwarzschild-ylciSs x 



2 J„2 

ds^ = -r2(l + f)dx+dx- + V(l - mdx+f + {dx-f] + r\dx^f + {dx^f] + ^ + dnj 

2 v^j 

G^,{X)dX^'dX'', (2.5) 



where / = 1 — dO,^ is the line element on and x^ corresponds to the light-cone 
coordinates on the boundary (r — )■ oo). In the last step of (2.5), we have introduced the 
compact notation of the target spacetime coordinates 

Xf" = {x^,x-,x^,x'^,r;n5) (2.6) 

with ^5 the set of five coordinates on 5^. 

The bosonic part of (2.2) ^^[X] is the Nambu-Goto action of a 2d world sheet embedded 
in the target spacetime, given by 

SbIX] = [dr [ dr^g, (2.7) 



2vra' 

where g is the determinant of the induced 2d metric 



dX^" dX 



V 



gai3 = Gf,^ {X) , (2.8) 

with a"{a = 0, 1) the world sheet parametrization. The fermionic part of (2.2), Sf[9], takes 
the form [23-25] 

SF[e] = ^ / d^cjiiv^g^^ep^D^e - ie^^^ep^ppO) (2.9) 
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where and are two 16 x 1 Majorana-Weyl spinors in the target spacetime and 

^X^' 



with Fa (a = 0, 1, 2, 3, 4) the 16 x 16 Gamma matrices satisfying 



(2.10) 



(2.11) 



and r] = diag(— 1, 1, 1, 1, 1). The higher powers of has been ignored and the bosonic coor- 
dinates in Sf\0] has been approximated by the classical solution. The Majorana condition 
relates 9 and 9 via 

9 = eCxC (2.12) 

with C the 4x4 charge conjugation matrix. The K-symmetry of the original fermionic action 
up the the order 6^ in the Schwarzschild AdS^ x background has been employed to reduce 
the fermionic degrees of freedom by half. We refer the interested readers to [26] for the 
detailed reduction for the present application. It is a straightforward generalization of the 
K-symmetry in the AdS^ x background [27-29]. 
The classical world sheet corresponds to 



= (0,r,0,x^r, 0,0, 0,0,0) 



and ^ = with the function x'^ given implicitly by 



dx'^ _ _l_ e 
dr r'^y/f' 



where 



{2-k)2TL 



(2.13) 



(2.14) 



(2.15) 



and the positive(negative) sign refers to the upper(lower) branch in FIG.l. The projection 




Figure 1. The curve represents vs r. 
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of the classical word sheet onto is a point. The induced world sheet metric 

ds^ = Qapda'^da^ = ^dr^ + ^^^dr'' (2.16) 

with a = (r, r), and 

= G^^^^- (2-17) 
and the world sheet scalar curvature reads 

^ = -rT^a-7r)' (2-18) 

and one can readily verify that 

j d^Gy/lR = 0. (2.19) 



The classical action reads 



S\X, 0] = Sb\X] = i — ^ VAL_tVi + (2.20) 

whose term of the Taylor expansion gives rise to (1.2) following the definition (1.1) 
(< H^^[C] >~< W[C\ >^ at large Nc). Because the metric of the classical world sheet (2.16) 
is space-like, its Nambu-Goto action and the quadratic action of the fluctuations around it 
are all imaginary. 

For the rest part of this section, it is convenient to introduce the funfbeins pertaining 
tangent space of the Schwarzschild-^dS's sector 

= r^{dX+ + dX-) 
E^ = -^(dX+ - dX-) 
E' = rdX\i = 1,2 

E^ = -^dX\ (2.21) 

which can be further decomposed into transverse and longitudinal components via a rotation 

£^ = E^ cosh, (i + E^ sinh p 
= E^ sinh 13 + E^ cosh f3 
S' = E' 

S"^ = e"^ cos a — E^ sin a 

8'^ = E"^ sin a + E^ cos a (2.22) 

with tanh /3 = ^/f and tana = e. 

The action underlying the partition function consists of the the quadratic terms of + 
6X, 9] in 6X and 9 with 

dXf'ia) = i6X+,5X-,5X\5X^,6X'^,6n5). (2.23) 
The linear terms vanish because of the classical equation of motion. 
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2.1 Bosonic fluctuations 



Among five fluctuation coordinates within Schwarzschild-AdS's sector, only three of them are 
physical and the rest two correspond to a world sheet diffeomophism. This can be made 
explicit be projecting 5X onto the funfbeins i5's, i.e. z"" = E'^6X^^{a = 0,1,2,3,4) and to 
retain only the transverse fluctuations and z^ and z'^. As the longitudinal fluctuations 
z^ and contribute only zero modes of the quadratic action, it is more convenient to work 
within the static gauge, 6X~ = SX'^ = 0, through a superposition of the transverse and 
longitudinal ones. Then G^u is independent of 5X and we have 

^ ,-,d6Xfd6X'' , 
= 5a/3 + G^.(X)^-^^-^. (2.24) 

Upon the parametrization 



u = bX^ = -jL-l^,v = 6X'^ = ^^^''\ ],w = 6X^ = -C (2.25) 
r y 2/ r r 

and (51^5 = 6s with s = 1, ...,5 with Xs the five tangent coordinates in the neighborhood of 
the classical point. Expanding the bosonic action (2.7) with -y/— 5 = ^J —grrgrr + Qrr to the 
quadratic order in 5X^ , we find that 

'^bosonic 7i 7 

(S<; + S^r, + Ss) (2.26) 



27ra' 



where 



s<: = lij dTdrv^ir^c' + rc'' + M^e) = zc>c (2.27) 

S^, = Uj drdr^gi-r^e - + g^^if + " M^f + - -^=^^^v) 

1 ~ /e 



= U{lfi)AA^] (2-28) 



and 



S5 = \ij drdr^g (5^-5/ + f'^) (2.29) 

where g'^'^ = ^,g^^ = F = ^ and F' = The diagonal masses in (2.27) and 

(2.28) read 

M,^ = M,= = s; + Mr7i = _^l_ = _M^ (2.30) 



and 



Notice the negative sign in front of the terms quadratic in ^. 
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2.2 Fermionic fluctuations 

To simplify the fermionic action (2.9), we take the representation 

= 7" X h (2.32) 

with 7" 4 X 4 gamma matrices satisfying 

{^-,^^} = 2t' (2.33) 
Upon a similarity transformation 9 = Vip with 

y = e-i/3ror3+|ar2r4 (2.34) 

we have 

V-\g'^p,D,)V = T^dr - -u;^r'^3T^) + T'dr + ^i==ror2r3 (2.35) 

where 

= -e^^Fg r = e^Ti (2.36) 
with the world sheet zweibines and spin connections given by 



^2 

.0 _ 'h 



^ V2r ' r\ f 

and 



(2.37) 



The fermionic action in terms of ip reads then 

Sf{0] = a /" d^a^i'lr^dr - loj^r'^sU) + T^dr + iTsU + ^i==ror2r3]V'. (2.39) 

Tra y 2 V 1 + 

The next step is to decompose the 16 x 1 into eight world sheet Majorana fermions. 
This is accomplished with the representation 

( -ia2 \ ( (Tz q\ 
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with ai{i = 1,2,3) Pauli matrices. It follows that 

\ asxhj y -aixhj 

r3r4=f-^"^''^ ° ] W3=f-"^"'^ ° ] (2.41) 

I -ia2 xh \ as X /4 i ^ ^ 



and 



where 



Sf[x\ = a / c^^c^^V^ {x+sA+x+s + X-sA-X-s) , (2.42) 



and Xits's are two component Majorana fermions. 

In summary, the explicit form of the fluctuation action reads 

A5[e, r?, C, 5, X] = ^(-^C^ + S^ + Ss) + iSplx] (2.44) 

with 5^, S^^, Ss and S'Fie'] given by (2.27), (2.28), (2.29) and (2.42). The partition function 
of the fluctuations is defined to be the path integral of 

Z = const. J [d(][dr]][d(][dS][dxy^^^^''^''^'^'''^ (2.45) 

and the effective action up to one-loop order reads 

SeS.=S[X,0]-ilnZ. (2.46) 

As was pointed out in [30], the classical world sheet underlying (1.2) is the saddle point of 
the Nambu-Goto action. Therefore the bosonic part of the path integral (2.45) is problematic. 
In particular, the path integral 

f -j^ilv)AJ^] 

'dC][dr]]e'23 = [d^][dr]]e V"?/ (2.47) 



with A^ri defined in (2.28) does not exist for a real and positive a' because of the negative 
eigenvalues of A^rj- In order to make the leading order term (1.2) stable and to make the 
one-loop correction to be computed below meaningful, we define the path integral (2.45) for 
a' > as an analytical continuation from an imaginary a' . In terms of the spectrum of A^ri, 
{ujn}, and a' = is with s > 0, the path integral (2.47) becomes 

J [dC][dri]e'^ = Yl (^27re±^f oc |det-5^5^| (2.48) 
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with the plus(minus) sign corresponding to a positive(negative) ojn- Though we are not 
able to come up with a physical interpretation of this procedure, we observe that the sign 
ambiguity associated to the continuation of the square root -y/s = \J —ia' from the imaginary 
ol back to a real and positive one may be removed by renormalizing (2.45) by its value at 
e = 0. This amounts to subtracting from (2.46) the one-loop effective action of two isolated 
light-like straight lines. Since the jet quenching parameter comes from the term of the 
Taylor expansion of (2.46), this renormalization is unnecessary for our purpose. Choosing the 
constant in front of the path integral of (2.45) appropriately, we end up with 

det^^+detM- ^ ^2.49) 

4^^det^^^deti(-V2) 

and the finite 't Hooft coupling correction to the jet quenching parameter comes from the 
term of the Taylor expansion of 

In^ , ^ 

z = - hm -— , (2.50) 

where L_ defines the domain of r, i.e. < r < L_. 

Because of the cancelation between bosonic and fermionic degrees of freedom, S'eg is at 
most logarithmically divergent. The coefficient of the divergence takes the invariant form 

y dVVKci + C2R) (2.51) 

It is shown in the appendix A that ci = as expected. While C2 7^ in the static gauge, 
this is not bothering because of (2.19) for the world sheet under consideration. We will not 
address the subtleties brought about by the world sheet boundary [31] but are content with an 
explicit demonstration that the one loop correction to the jet quenching coefficient is indeed 
divergence free in the next section. 

3 The formula for k 

In this section, we shall transform the effective action (2.46) into the form suitable for the 
numerical calculations described in the subsequent section. In the course of the reduction, we 
shall demonstrate explicitly that the correction to the jet quenching parameter is free from 
both ultraviolet and infrared divergence. 

It is convenient to work with the new coordinates 



where the section X2 > 0(the upper branch of FIG.l) of the world sheet is mapped to < 
p < 00 and the section X2 < 0(the lower branch of FIG.l) to —00 < p < 0. We have 

pc.±^ (3.2) 



- 10 - 



for large r and 



1 

^ 2 



for r near the horizon. So the two branches of r to /3 mapping join smoothly at 
terms of the coordinates {t, p), the world sheet metric takes the conformal form 



with 



2r2- 



The corresponding zweibines and spin connections read 
and 



01 



r V 2 dp 



The functional operators underlying the determinants in (2.46) becomes 



where 



■Au 



1 + 



92 a2 



■Au 



9^ + 9^)+^'^ 



_4f:p-</'A 



-20 



dt"^ dp^ 
( d i d(j) \ d 



+ i(T2 Vl + =p ecJs. 



The partition function becomes then 



det^^+det^^- 



det 2 ^^^det 2 ^^det 2 ^„ 
Upon a Fourier transformation in t, we find that 



with 



det 2 (w ) det 2 (w ) det 2 (i^ ) 

and z = ^ ^j^lnZ(w), where ^...(f^) is obtained from A,,, of (3.9) by replacing ^ with 
ExpHcitly, we have 

Ai:^,{^) = e-^^D^r,{oo) 
Auiuj) = -e-^'t'D^iuj) 
A±{uj) -- 



(3.12) 



-lUJ. 



where the functional operators 

^o(^^) = 



dp' 



D(-r,{uj) 



d 



-cji— + iwfJs + e''(iv 1 + £^^2 =F ecis) 
dp 



with 



(3.13) 



(3.14) 



dp' 



+ a;2 + e^'f'M^ 



Taking the square of D±{(jj), we find 



d^ 
dp^ 

-V±V±{oj)V^^ 



=F 2iea;e'^ - e^'*' + (Vl + eVs T iea2)ef 



(3.15) 



dp 



(3.16) 



where V± is the 2x2 matrix that diagonalizes the matrix Vl + e^cs =F ^£'^2 and 

ct/o"^ dp 

_ / Dsiuj) T 2i€uje'l' \ 

~ I D'r^ioj) ^ 2ieu}e't' ) 



(3.17) 



with 



2(/) 



(3.18) 
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Replacing the p in D-i{uj) by —p, we can get D'^{(jj) and therefore detL'3(aj) = detD3(a;). 

The scaling factors e"^''^, e""^, e~2'^ and ea*^ of (3.13) cancel in the ratio of determinants 
of (3.12). Consequently 



In the limit 



oo, 



detP+(a;)detP_ (w) 



det 2 D^^{oj)det 2 Di (u;)det 2 Dq {< 



Loo27r^ > 



(3.19) 



(3.20) 



with L_ 



7==^L_ the length of the domain of t. It follows that 

T 



z 



dw In Z{oj). 



Since Z{oj) is an even function of e, we have the Taylor expansion 

In Z{uj) = a{uj) + b{uj)e^ + O(e^) 

and 



T 
'2 



T 

doja{uj) + e^ — 



-a{Lo)-h{uj) +0(e^). 



(3.21) 



(3.22) 



(3.23) 



It is the second integral that contributes to the correction to the jet quenching parameter. 
Hence, the strong coupling expansion of the jet quenching parameter can be written as 



where 



80F 



r2(i) 



doj 



(3.24) 
(3.25) 



3.1 The formulae of a{uj) 

The function a{ijj) corresponds to lnZ(a;) at e = and we have 

det^L>3(a;) 



a{uj) = In 



detl^i (tj)det 2 D2 (a;)det 2 1)0(0;) 



(3.26) 



a{uj) can be factorized as a product of determinant ratios, which can be evaluated by the 
Gelfand-Yaglom method [32] [33]. We give a brief introduction here. The examples of its 
application to other Wilson loops can be found in [34-37]. 
Consider two functional operators 



Hn 



--^ + Vn{x), 



(3.27) 
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with n = 1,2, defined in the domain a < x <h under the Dirichlet boundary condition, the 
determinant ratio 

deti/2 f2{h\o) 



detFi h{h\a) 
where fn{x\a) is the solution of the homogeneous equation 



(3.28) 



Hnfn = 0, (3.29) 

subject to the conditions /n(0|a) = and /^(0|a) = 1. In terms of a pair of hnearly indepen- 
dent solutions of (3.29), {rjniin) 

where the Wronskian ^^[r^njCn] is ^-independent. 

For our problem, we have Hi = Dq{uj), Di{uj) or D2{uj) and H2 = D^{u)) under the limit 
a — )• — 00 and 6 — )• 00. Let us define by u~a^ and u~a'^ the two linearly independent solutions 
of the equations -DqV = 0(a = 0, 1, 2, 3) subject to the boundary conditions 



e '^P p ^ 00 



We have Co{uj) = 1 and the rest of Ca{uj)^s satisfy 

= -2^^C„(w), (3.32) 

It follows from (3.26), (3.28) and (3.30) that 

0(0;) = In = 41nC3(a;) - lnC7i(a;) - - lnC2{uj) (3.33) 

CliL0)CiiL0) 2 

The asymptotic forms of Ci{lj), C2(ct;) and 6*3(0;) for large lo and small lo can be determined 
analytically. We find 

: n f ^ 0-61802 , ^ , , 5TH\) 6.1802 0.92704 
lnCi(w)~— lnC2 6j lnC3(w)~-^~ , 

(3.34) 

as — 7- cxo and 

n ( ^ 0.927 ^ , , 225r^ (|) 26.0736 1 

as w — 7- 0. The details of the derivation are shown in the appendices B and C. It follows from 
(3.33) that the integration over uj is convergent. 
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3.2 The formulae of b{uj) 

With a nonzero e, the Gelfand-Yaglom method becomes comphcated as the differential equa- 
tions of and 7] are coupled, but we could employ the perturbation theory to calculate b{uj). 
From (3.19), we have 

1 15 

In Z{u}) = 21ndetP±(a;) + -\ndetD^,^{uj) - -\ndetDi{uj) - -lndetDo(w) (3.36) 

where 'D±{uj) is given by (3.17) and D^fj{tj), Di{uj), Dq{oj) are given by (3.14). For conve- 
nience, we denote 

1 1 ~ 

Z^rii^^) = --IndetD^rii^) = - ^ l^det[£'g,,(a;) + 6D^ri{^)] 

Z±{uj) = 21ndetP±(w) = 21ndet[P±(a;) + 5V±{uj)] (3.37) 

where 

^.H=f''^^"^ ] SV,i.)=(^'''^'' ' ,] (3.38) 

with D3{u}) and -03(0;) given by (3.18) and D2{oj) given by (3.15). The next step is to expand 
Z^^{uj) and Z±{uj) to the 2nd order in 6D^^{uj) and 6'D±{uj). The 1st order terms do not 
contribute because ti6D^^{Lj) = ti6'D±{Lj) = 0. Consider first the Z^j^{uj) 

Z^r^^uj) = —- In det 1)^^(0;) 

= ~ lndet[L>g^(a;) + ln(l + SD^^I{uj)6D^^{u}))] 



-^lndet55,(a;) - ^TrD^;;{u)6D^^{u) + ^TtD^^';{u)6D^^{u;)D^^I{uj)6D^^{u) 
Z^r,{uj)+6Z^^{oj). (3.39) 



Likewise, 



Z±{uj) = Z±{lo) + 6Z±{u}). (3.40) 
One can readily verify that at e 7^ 

In Z{uj) = a{u) + 5%(a;) + 6Z±(u}). (3.41) 
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where a{u)) has been discussed in the previous section. It then foUows from (3.22) that 



(3.42) 



where 



-D^\u) \ / 
D^\io))\- 



Aieuoe'i' \ ( -D^^iuj) 



Aieue'^ 



dp / dp'gUP - P)e'^^''^92.ip - />')e*(^) 



where the Green functions giuiip — p'), g2u}{p — p') satisfy 



dp' 
£_ 
dp'^ 



(-— +^' + e'''M^)gUp - P') = Hp - p') 



+ a;2 + e^'^Ml)g2Up " P) = Kp " P)- 



In terms of the solutions defined in (3.31), we find 



giu{p- p') 



g2Lo{p- p) 



1 



2ujCi{u}] 



1 



4'J(p')^S2(/') 



(3.44) 



(3.45) 



(3.46) 



2ojC2{oo) 

for p' < p, where the dependence on lo is indicated exphcitly. In what follows, we shall 
suppress the superscripts and denote u^alip) by Uaui{p) for a = 1,2. It follows from the 
symmetry of their differential equations that UaMp) = Uauj{—p)- 
Substituting (3.45), (3.46) into (3.43), we have 



4e^ 



dpe'i'^P^u^^{-p)u2^{-p) / dp'e<^(^')ui,(p')^^2c.(p)-(3.47) 



Ci{uo)C2{lo) 

The process to calculate 5Z±[uj) is similar to the 5Z^r}{'^) part, here we just show the result 



5Z±(w) 



4e^ 



dpe'^^'^vlip) / dp'e'^^f'^uUp'). 



(3.48) 



where U3^{p) and V3aj{p) stand for ul^\p) and u'^\p) defined in (3.31). Applying (3. 47), (3. 48) 
into (3.42), we find 



+ 



Ci{oo)C2{oj) 
4 



C|(u;) 



dpe'^^P^ui^{-p)u2^{-p) / dp'e^^p'^uUp')u2M 

D J — OO 

dpe^^^^Up) r dp'e^^'^'^Up ). (3.49) 



(3.43) 
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To prove the convergence of the integral of b{uj) over cj, let us examine the large lo behavior 
first. The solutions of -01^2,3 (f^)""!, 2,3 = may be approximated by ni^2,3(p) — e'^'' and 
Ci,2,3 = 1 + izj)- Introducing the new variables X = \{p + p') and x = p — p' , We have 
then 

dpe'^^P\U-p)u2u.{-p) r dp'e'^^P'^uUpMp) 



Ci{oj)C2{oj) 

r dX r dxe<^(^+f )+<^(^-f = -1 r dXe^'^^^^^ + o [1] (3.50) 



and similarly 



2 



dpe^^'^^vlip) r dp'e^'^P'^Up) 



Consequently, the leading terms of the two integral in (3.49) cancel and b{uj) vanishes faster 
than O as a; —7- 00. 

The integral of b{uj) is also convergent at the lower limit because 

lim b(uj) ~ 2. (3.52) 

the details are shown in appendix D. 
4 The numerical calculations 

In this section, we present the numerical calculations of a{uj) and b{uj). 
Consider the equation 

dhi 
dp' 

where V refers to the potential of the Schrodinger like equation DiU = 0, (i = 1,2,3). To 
avoid the exponentially growing behavior, we factor out e^^ and set u{p) = r]{p)e^P, then the 
equation (4.1) becomes 

dp"^ dp 

which is suitable for the numerical calculations, then coefficients, Ci{uj), C2(w), C^{uj) to- 
gether with the solutions r]i{p) can be obtained from the equation (4.2). This method has 
been employed recently in Ref [35] and [36]. 



2 {uj' + V)u = ^. (4.1) 



2 + 2w-j^ - l/r? = 0. (4.2) 
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We start from the negative value of p with \p\ » 1, where these Schrodinger Hke equa- 
tions take the asymptotic forms 

(f 

Diui ~ {--r^ + uP')ui = 
D2U2 ~ ( — + ^ + —^)u2 = 



2 

D3t/3 ~ (-— + a;2 + — )n3 = (4.3) 



At p = —K {K >> 1 is a large cutoff, and K — )■ 00). We can find the solutions of (4.3) 

ui{-K) = e-^^ 



U2{-K) = (1+ + 



oK 



u,{-K) = (1 + -L)e-^, (4.4) 



It follows from ri{p) = e ^Pu{p) that 

mi-K) = 1, v'ii-K) = o 



which serve as the initial conditions for the numerical solutions. Then we run the forth order 
Runge-Kutta algorithm all the way to p = K. There, we can find the approximate solutions 
as well (See appendix C for details), from which Ci{uj), C2{oj) and 6*3(0;) are extracted with 
the formula 

Ci(u;) = r?i(i^) + ^^(^) 



C2{U 



2uj 

ri2{K)N + [i^{K)+u:r,2{K)]Q 
2uj 



where 



C3(a;) = m{K) + ^ (4.6) 



3 6 6 

N = 0J+ — + ——: + 



Q = l + 4^ + ^. (4.7) 

The numerically generated Ci{uj), C2(w) and 6*3(0;) for large uj and small u read 

, „ , , 0.61799 , ^ , , 6.1799 , ^ , , 0.92703 ,^ ^, 

lnCi(a;) = lnC2{u}) = In 63(0;) = (4.8) 

UJ UJ OJ 
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0.927 



26.0971 



0.5001 



(4.9) 



Comparing the above behaviors with (3.34) and (3.35), we find that the numerical results are 
in agreement with the analytical ones. Coming back to (3.25), owning to the evenness in ui, 
we have 



did 



]^a{uj) - b{uj) 



16V^ 







°° r 1 

dhj —a{(jj) — b{u}) 



(4.10) 



Having obtained the numerical values of Ci(a;), 6*2(1^), C3(tj), the integration of ^a(a;) is 
carried out in the following steps. The Simpson rule yields 



.14 ^ 

/ du-a{u) = 1.548. 
Jo. 001 2 



(4.11) 



Considering the small u behavior of Ca from (3.35), we can fit the expression of ha{oj) 



-a{oj) ~ -2.155 - 2.249 ln(w) 



as cj — )• 0. Thus 



0.001 



duj—a{uj) 



0.001 



da;[-2.155 - 2.249 ln(a;)] = 0.016. 



Likewise, the expression of ^a{uj) when uj is large can be fitted to 

1 , , 0.002 0.646 
a(w) ~ — ^ h 



then 



14 



1 , , /■'^ , ,0.002 0.646, 
du;-a{Lo) ~ / dio[ — ^ — I 3— J 



0.002 



Combining (4.11), (4.13) and (4.15), we have 



00 -j^ 

dbj—a(oj) = 1.566. 



The function b[uj) can be divided into two parts 
4 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



+ 



Ci{oo)C2{oj) 
4 



dpe'^^p\i^{-p)u2^{-p) / dp'e^^p'^uUp')u2M 



dpe^^'^vUp) / dp'e^^p\l{p') 



(4.17) 



where 61(0;) and 62 (i^) are the 1st and the 2nd double integrals on RHS. In terms of the 
numerical values of Ci(c;j), 6*2(0;) and 6*3(0;) and the corresponding solutions i]iujip) at each 
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5 1 1 5 2 2 5 (j^ 3,0 



Figure 2. The curve of vs uj. 



oj, the numerical values of 61 (w) and 62(0;) can be obtained. Combining bi{u}) and b2{io), the 
curve of 6(0;) is shown in FIG. 2. Notice that the integral of b{oj) over u is convergent since 
b{Lo) falls off faster than uj~^ towards the upper limit and approaches a finite limit at lower 
limit, consistent with the analytical results of last section. When u is not large, we compute 
the integral by means of the trapezoid method and obtain 



/ dujb{uj) ~ 2.283. 
Jo 



For large u, b{uj) can be fitted to 



0.21 2.329 



and the integral is carried out analytically. We find 



^, , f°° ^ ,0.21 2.329, 
du}b{u) ~ / du[ — ^ H ^ 

3 ^3 ^ 



0.197 



Combining (4.18) and (4.20), we have 

dwb{uj) = 2.283 + 0.197 = 2.48. 

From (4.16) and (4.21), we find 

/ duj[-a{uj) - b{uj)] = 1.566 - 2.48 = -0.914. 

Substituting (4.22) into (4.10), we end up with the numerical value of n 

16V^ 



00 -j^ 
dijj\—aiijj) 



b{u:)] 



-1.97. 



The physical significance of k will be discussed in the next section. 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



(4.23) 
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5 Concluding remarks 



The jet quenching is an important signal of the quark-gluon plasma created via heavy ion 
collisions. At the energy scale of RHIC, the QCD running coupling constant is not weak, 
the perturbative analysis may not be appropriate, nor the lattice is available. In terms of 
the QCD running coupling constant gqcB and oqcd = 4^° , the 't Hooft coupling Aqcd = 
47rA''caQCD = 127raQCD and is strong when aqcD = This motivates the application 

of AdS/CFT to the A/" = 4 super Yang-Mills plasma as an important reference. To the 
leading order of large Nc and strong 't Hooft coupling A, the jet quenching parameter in the 
super Yang-Mills plasma has been calculated. The sub-leading term in the strong coupling 

3 

expansion, which dominates the 0(A~2) term in the literature, is provided in this paper. 

To figure out the leading order correction, we have expanded the Nambu-Goto action 
around the word sheet underlying the order jet quenching in the quadratic order in the 
fluctuation coordinates and carried out. As the leading order world sheet is a saddle point of 
the Nambu-Goto action, we defined the one loop effective action as an analytical continuation 
from a path integral with an imaginary string tension, from which the correction is extracted. 
The resultant effective action is expressed in terms of ratios of determinants and is shown 
explicitly to be divergence free. The reason for the absence of UV divergence can be attributed 
to the the metric (2.5) solves the lOd supergravity equation of motion [25]. The numerical 
computation of the determinant ratios yields the result (1.4). The negative sign of k is 
consistent a monotonic behavior from strong coupling to weak coupling. For asYAl = 1/2 
and Nc = 3, we have A = 18.8, which gives rise to 32% reduction of qsym from the leading 
order amount. 

The jet quenching parameter has also been calculated within the AdS/CFT correspon- 
dence via the dragging force with a slightly different outcome to the leading order. As this 
calculation also involves the minimum area under the Nambu-Goto action, it remains to ex- 
amine the order A~2 along this line. This discrepancy from the experimental result may also 
stem from other difference between QCD and the super Yang-Mills at large Nc and large 't 
Hooft coupling, such as finite Nc and infrared cutoff of the former. We hope to report our 
progress in this regard. 
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A The UV divergence of the bosonic and fermionic fluctuations 

In this appendix, we shall prove the cancelation of the UV divergence caused by the mass 
terms of the fluctuation action (including the cross term of (2.28)), which contributes to ci 
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Define 



of (2.51). The world sheet curvature does not contribute in this case and we may work in an 

cartesian frame in the neighborhood of each point, do = (ro,To). To make ah path integrals 

meaningful, we shall start with an imaginary a' in accordance with the discussion following 
(2.49). 

1. The C mode 

Sc = y drdr^girC" + ^ C'' + M^C^) (A.l) 

1 

r — ro = xi 

1 , 

r — To = _ — X2 (A. 2) 

ds^ = dx\ + dx\ 

j dT j dr^ = J d^x (A.3) 

Sc = ^ASc (A.4) 

Co 

= ^ y d^x{d^CdaC + Mie) (A.5) 



then 



and 



where 



Making a Fourier expansion 



ipr 



where is the area covered by the cartesian frame at (Tq, then 



V 

g = c-p (A.6) 



^^c = ^ T.^f + M'^KpCp = Y^{f + mDqc,^ (A.7) 



p 



with extending to half momentum space. 
On writing 

j \\_dC*;dC^^^^^ = const.e-^^S (A.8) 

p 

where 

Ar, = i ln(p- + Ml) = \ Inf + hl^ 1 (A.9) 
p p p 
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and the constant does not contribute to q. So the mass induced logarithmic divergence 

(A.IO) 



(UV)c = -M2 —Mf /"^ ^ ^ —M} In A 

V >^ ') ( p2 47r ^ y p 47r ^ 



p<A ' 



2. The and r] modes 



8\/2, 



e r 



where 

Making a Fourier transformation 



ipr 



then 



On writing 



we find 



Y{ dCpdipdr]*;dri^'2^^^'i-^ = const. e"^'^^'' 



„2 I /Vf2 _ 4iep2 



Then the mass induced logarithmic divergence 

p<A 



(A.ll) 



id2v) (A.12) 



(A.13) 



p 

J^[{f + M|)e;^ep + (p^ + M2)r?tr?p^+ -^i=p2(e;^r?p - r?;^^)] (A.14) 



(A.15) 



(A.16) 



(UV),, = i(M| + Ml + E ^ - ^(^1 + + 1- A (A.17) 
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It follows that the coefficient of the total logarithmic divergence due to the mass terms of the 
bosonic fluctuations 

= 11 d^<^V^{8 + R) = ^j ^'^VI (A.18) 

3. Fermionic modes 

The operator (2.42) can be related to a massive Dirac operator of in a 2D curved space 
and a U{1) gauge potential, i.e. 

^± = ZCJ2 (7"V^ + m) , (A.19) 

where the mass m = 1 and the covariant derivative 

Va = ^-iAa + lb]ujt (A.20) 

with the 2D gamma matrices 7*^ = cri, 7-*^ = —(73 and the gauge potential Ar = =bi-7=5e°, 

Ar = 0. As is in a flat 2D space, the gauge potential does not contribute to the UV divergence. 
Therefore, the UV divergence of the fermionic effective action in (2.49) 

= 2(lndet^+ + lndet^_) (A.21) 

is the same as 

r'^ = 4 In det^ = 2 In det^^ (A.22) 
where A is obtained from A± with the covariant derivative replaced by 

Va = ^ + |[7,76K'. (A.23) 

We have 

^2 = v"V„ -1-^R. (A.24) 
It follows that the coefficient of the logarithmic divergence is given by 

J dVV5(4 + const x i?) = 4 J dV^I, (A.25) 
which cancel the bosonic contribution. 
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B The large oj behavior 



In this appendix, we present the details of the derivation of Ci{ijj)^C2{ijj)^C^{ijj) at a large 
and positive oj. 



We start with the 2nd order ordinary differential equations 

dp 
^ dp 

and 



(--^ + ^2^e2'^Mf)/i = 0, (B.l) 



2 + ^2 + e^^Ml)f2 = 0. (B.2) 



--^ + ^2 + e2^-e^j^)/3 = 0. (B.3) 
ap^ dp 



Next, we find the WKB approximation large uj 



and 



j,^^±Fdp'^^^-e^t (B.5) 

. By choosing appropriate sign and the integration constant pQ in each case, we may find the 
large oj a 
we have 



large uo approximations of the specific solutions Uq ^ and Ua* defined in (3.31). In particular. 



and 

Thus, we find that 



^(l)(^)^,-P+^/_^.^p'e^^M^^ (B.6) 

nf (p)-e-^+^/-o^'^'''^''^= (B.7) 



lnC^{u)^^ dp^^Ml = \^ 



1 r .__2<*,.2 5r^(i) 



1 r°° 

lnC2(a;) ^ — / dpe^-^M 



lnC3(a;) - — T ^^(e^* - e*^) = (B.9) 

as cj — )• oo. The final results in terms of the Gamma functions are obtained by transforming 
the integral with respect to p back to that with respect to r according to (3.2). 
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C The small ui behavior 



To estimate of the small uj behavior Ci, C2 and C3, we introduce three domains of p, Ir l '■ 
\p\ « 1 and II : \p\ « 1/uj with the subscripts R and L referring to the right side {p > 0) 
and the left side (p < 0) of the origin. For a; << 1, the intersects II n II and // H Ir are 
nonempty and serves as buffers where the analytically approximate solutions within different 
domains can be matched. In what follows, we shall denote the set of functional operators 
involved hy D = {Di, D2, D3, D'^} and their kernels by u = {ui, M2, ns, 1*3}. Though the 
subset of {Di, D2, -D3} is sufficient to determine Ci{uj), C2{uj) and 6*3(0;), we include D'^ in 
the early part of this appendix since v^^^lp) of (3.31) can be obtained from the solution of 
Dgtig = subject to the same boundary conditions of ui, U2 and upon replacing p by -p. 

Within the domains II and Ir, all powers 1/p in D that fall off faster than l/p^ can 
be ignored and the equations D(j) = become analytically soluble in terms of elementary 
functions. We have 

2 

D2^-^+u' + ^ (C.l) 



for both II and Ir, and 



for II and 



dp2 A dp^ p2 ^ I 

for Ir. In the domain //, on the other hand, w can be approximated by zero and -D's become 

ap^ 

ap^ dp 
(ip2 dp 

In the overlapping regions II H // and // n Ir where both high powers of ^ and cj can be 
dropped, we end up with 

^ " dp2 
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for both Ii n // and II fl Ir, and 

dp'^ ' p'^ ^ dp 



within II D II and 



within II n Ir. 

In what follows, we shall start with the approximate solutions in II with the asymptotic 
form e''^P as p — )• — oo, matching different approximations through // and Ir and obtain the 
leading small oj behavior of the coefficients of the asymptotic form e'^P as p — > oo. 

The solutions of Du = in with D's given by the approximation (C.l) and (C.2) 
subject to the condition that u — )• e'^'' as p — )■ — oo read 



I 2ujp 
U2 — \j Ks{—ojp) 



U3 - \l -'^^Ksi-ojp) 
V vr 2 

u's ^ e^P, (C.8) 

where Ky{z) is the modified Bessel function. In the left overlapping region I^ n // where 
—ujp « 1, the solutions (C.8) can be approximated by 

Ml ~ 1 

3 

U2 ~ 



1 

M3 ~ 

cup 

n'3 ^ 1, (C.9) 

matching the power law behavior dictated by (C.5) and (C.6). It follows from (C.5) and (C.7) 
in the right overlapping region that 

Ml ~ + Bip 

1 ( A2 . 



ufi \ p2 

M3 ~ -(^3 + S3P) 

+ (C.lO) 

LjJ p 

where the coefficients A's and B's are determined by (C.9) and the solutions covering the 
entire domain II, with L)'s given by (C.4) and are independent of w. Coming to I/j, the 
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solutions take the approximate forms 

ui ~ 2(Cisinha;p + (7ie-^'') 



U2 ^ v27rcj/j[C'2/5 (wp) + C2K5{ujp)] 

2 2 

U3 ^ 2(C3 sinhwp + Cse^'^P) 

U3 ~ y^2TT0jp[C3l3{ujp) + CgKlCwp)], (C.ll) 

following (C.l) and (C.3). All these solutions grows exponentially as p — t- 00. The reason 
why U3{p) and u'^{p) share the same coefficient of exponentially growing follows from the fact 
that u'^{—p) is the kernel of D3 and its Wronskian with U3{p) is independent of p. Matching 
the approximate solutions of (C.IO) and (C.ll), we find that 

61 ~ — 

^ 15^2 



ft |L (C^12) 



and 

(7iH = 0(l), (72(a;) = 0(1), C3 = 0{u:'^), C^(a;) = 0(a;) (C.13) 
for u « 1. 

The coefficients i?'s can be figured out analytically. It is more convenient to work with 

the old variable r defined by the transformation (3.2). A large r corresponds to a large 
magnitude of p according to 

pc^±^ (C.14) 
with positive(negative) sign for p > 0(p < 0). Near the origin p = 0, we have 



as /? — 7- 0^. In terms of r-coordinates, we have 

and the equations Diui = and D2U2 = can be reduced to hypergeometric equations with 
the solutions that matches the powers of (C.9) 
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for p < 0, where z = ^ and the asymptotic form (C.14) has been used for matching. It 
fohows from (C.15) and the formula 

F{a,b;c;z) = ^^f}I^-^—^F{a,b;a + b-c+l;l-z) 
r(c — a)T[c — b) 

+ ^^^Y^^W^^^ ~ -a,c-b;c-a-b+l;l-z) (C.18) 

that the analytic continuation of (C.17) to p > reads 

u, = -F('-l'--,z)+'T^a)z~l 
V2 4'4' J 2V2^ V4y 

3 1 / 3 9 \ 15 fl\ 3^/1 1 1 \ , 

U2 = —.z-2F{i,-;-;z] + —^T^ -] z"^ F {--,--;--■ z) . (C.19) 



2t^2 V '4' 4' J 8V2^oj^ 
Matching (C.19) to (C.IO), we find that 



which yields the small lo behavior 



r2(i) 225r2(i) ^ ^ 



following from (C.12). 

As to C3, we notice that D3 of (C.4) factorizes as 



A + e-^ ) = (C.23) 



The solution of the 1st order equation 

also solves the 2nd order equation D^tp = 0. In terms of the r-coordinate, eq.(C.23) takes the 
explicit form 

7-f ± = (C.24) 

ar r J 

with the upper(lower) sign for negative(positive) p. Its solution 



matches the power behavior of U3 in the left overlapping region, (C.9) and joins smoothly at 
p = 0. Therefore we identify ^3 = ip. Matching U3 to the behavior (C.IO), we obtain i33 = 1 
and then 

ft ^ (C.26) 

as w — )• 0. 
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D The limt lim^^o K^) 

As w — 7- 0, we may find introduce po such that po >> 1 but \ujpq\ « 1 and divide the domain 
of the double integrals for b{uj) into the following six parts: 

I ■ -po < p < Po, -po < p' < Po, p' < p; 
II : Po < p < p; 

III : p> Po, -po < p < Po; 

IV : p> Po, p < -po; 

V : -po < p < Po, p < -po; 

VI:p'<p< -PO. (D.l) 

Correspondingly, we write 

b{uj) = bi{uj) + 677(0;) + 67//(w) + biv{uj) + bv{uj) + bvi{uj). (D.2) 

If — Po < P < Po, the solutions ui, U2, U3 and ^3 of the integrand may be approximated by the 
kernels of (C.4). Otherwise, the asymptotic forms (C.8) or (C.ll) are substituted for p < —po 
or p > Po- It follows from the approximate solutions in the appendix C that 

po rpo 

dpe^^P^uU-p)u2U-p) / dp'e'i'^P ^uiUp>2Up') = 

po fPO 

dpe't'^'^vlip) / dp'e^^P KlJp') = 0{oj~') 
-po J~po 

00 rpo 

dpe^^P^iU-pM-p) / dp'e^^P'^u,Up)u2Up) = 
Po -^-po 

dpe^^'^^vUp) / dp'e^^'^'KUp') = 0{u-' Inu) 

po J~po 

dpe't'^'^^uiU-pM-p) / dp'e^^f' ^mUp'Mp ) = 0{uj-^) 
po J -00 

00 y Po 

dpe^^^'^Up) / dp'e^^'^'^ulip') = 0{u-' Inoj) 

po J —00 

po r-po 

dpe'i'^P^uU-pM-p) / dp'e'^^P ^uUp>2Up') = 
-po J -00 

00 y po 

dpe'^^P^vlip) / dp'e'^^P )^L(p') = 0{u-') 
po J —00 



(D.3) 



Upon multiplication with C's, we find that 



lim bj{uj) = lim bjjj{uj) = lim biv{^) = li™ bv{uj) = 0, (D.4) 

0J—>-0 uj—^O uj—>-0 uj—^O 
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because of (3.35). As to bjj{uj) and bvii(jj), we have 

rp 



pa -J po 

Po , fP 



oo 

°° dr f'^ dv 

—e-^K,{x) / ^[Ci(a;)smhy + (7i(^^)e-^][C2(^^)/5(y) + C2Mi^5(y)] 



dpe'^('')ni,(-p)t/2.(-p) / dp'e^^P'\i^{p')u2^{p') 



p 



PO J PQ 

4 r ^[Cs{u)sinhy + Cs{u)e-y]' dpe^^^^Up) j' dp' e'^^^'^uUp') 



= 4/ dxKlix) dy[C-i{u)I-s{y) + cf\u)K.{y)]\ (D.5) 
It follows from (C.13) that 

lim bn{u:) = -16 r ^e-^Ks {x) {y) + 16 H — e'^- ^ sinh^ y (D.6) 

Jo ^ Jo Vy ^ Jo X Jo y 



uj^O 

and 



/■°° dx dv 1"°° 

lim bv I (uj) = -16 ^e-'^K'Jx) —hAy) + lQ dxKsix) dyldy). (D.7) 
Jo y/x ^ Jo y/y ^ Jo 2 _/o 2 



Therefore 



lim b{u}) 

LO^O 
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